Quantum analysis of shot noise suppression in a series of tunnel barriers 
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We report the results of an analysis, based on a straightforward quantum-mechanical model, of 
shot noise suppression in a structure containing cascaded tunneling barriers. Our results exhibit 
a behavior that is in sharp contrast with existing semiclassical models for this particular type of 
structure, which predict a limit of 1/3 for the Fano factor as the number of barriers is increased. 
The origin of this discrepancy is investigated and attributed to the presence of localization on the 
length scale of the mean free path, as a consequence of the strictly 1-dimensional nature of disorder, 
which does not create mode mixing, while no localization appears in common semiclassical models. 
We expect localization to be indeed present in practical situations with prevalent 1-D disorder, and 
the existing experimental evidence appears to be consistent with such a prediction. 
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In the study of low-dimensional devices, suppression 
of shot noise with respect to the value predicted (for the 
case of a Poissonian noise process) by Schottky's theo- 
rem has represented one of the most active fields of in- 
vestigation in the last two decades. Such a suppression 
phenomenon, described by means of the Fano factor, i.e. 
the ratio of the actual shot noise power spectral density 
to the full value 2q\I\ (where q is the value of the ele- 
mentary charge and / is the average value of the current 
fiowing through the device), has been predicted and ob- 
served in many different mesoscopic structures, and is 
the result of the presence of correlations between charge 
carriers, which reduce current fluctuations in the device. 
From the theoretical point of view, shot noise suppression 
has been investigated both with quantum mechanical and 
with semiclassical approaches. Such activities have led to 
the discovery of "universal" values for the Fano factor in 
specific structures; in particular, for disordered conduc- 
tors a universal suppression factor of 1/3 has been found, 
both with random matrix theory [l|, Q and with a semi- 
classical approach This result has received further 
confirmation from numerical simulations [3, 0] and from 
experimental evidence 6|| ■ 

A remarkable addition to these results is a derivation 
by de Jong and Beenakker 7, 's'l , who demonstrated that 
a Fano factor of 1/3 is obtained, within a semiclassical 
model based on the Boltzmann-Langevin equation, also 
for a series of barriers. A formulation relying on an equiv- 
alent semiclassical circuit model (similar to that used in 
for a series of chaotic cavities) leads to the same result. 

However, the only existing experimental data p^ . ob- 
tained for a GaAs/AlGaAs superlattice, are not in agree- 
ment with these semiclassical conclusions: they exhibit 
a Fano factor that depends strongly on barrier transmis- 
sion, in the limit of vanishing applied electric field. 

Prompted by this discrepancy, we have performed a 



quantum calculation of shot noise suppression for a struc- 
ture consisting of cascaded barriers. While for a conduc- 
tor with 2-dimensional or 3-dimensional disorder quan- 
tum simulations recover 0, [E| exactly the same 1/3 sup- 
pression predicted by random matrix theory and by semi- 
classical models, this is not the case in the presence of 
1-dimensional disorder, i.e. of randomly spaced cascaded 
barriers, regardless of the dimensionality of the conduc- 
tor. Our quantum calculations, although performed for a 
model system, do, instead, exhibit a behavior consistent 
with the experimental data. 

In this communication, we focus on the reasons for the 
discrepancies between the semiclassical and the quantum 
approach, and address the issue of what model best repre- 
sents the situation of a practical experiment. The struc- 
ture we have considered (a series of cascaded barriers, 
sketched in the inset of Fig. [1] and corresponding to that 
of Refs. 0j[ll) can be studied very straightforwardly from 
a numerical point of view. For each fixed value of the lon- 
gitudinal coordinate x, the potential profile is constant 
for Q < y <W {y \s the transverse coordinate, and W is 
the width of the structure: 8 /im in these calculations), 
while we assume a hard-wall confinement for y = Q and 
y — W . Thus, the orthonormal set of transverse modes 
is the same all over the structure, and the tunnel barri- 
ers do not introduce any mode-mixing. Since the modes 
are uncoupled, the numerical analysis can be split into a 
set of one-dimensional problems, one for each considered 
mode. For the generic n-th mode, the S-matrices of an 
interbarrier region of length L and of a tunnel barrier of 
height U and length I are: 
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where (at the Fermi energy Ep) 

A-2 _ 1.2 

Tn = exp(ifc„L) , pB„ = i -777-^ sin(/cB„0 tb„ 
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FIG. 1: Fano factor for a series of identical barriers, as a func- 
tion of the number of barriers for the case of realistic barriers 
(solid lines), and for model barriers with a transparency in- 
dependent of the longitudinal electron wave vector (dashed 
lines). Inset: sketch of the analyzed structure. 



with: k„ 



i-2 _ i.2 



kp — ^/2rnE^/h , kT„ = mr/W and ku = V^rnU /H . 
For each mode, the scattering matrices of adjacent slices 
are recursively composed to find the overall S-matrix and 
in particular one of its elements, the transmission i„ of 
the n-th mode through the device. The conductance and 
the shot noise powe r sp ectral density are then computed 
using the relations [ll|, [iJl 
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r„ (1 - T„) , (3) 



where the sums are performed over all the N modes prop- 
agating in the interbarrier regions, r„ = |t„p and V is 
the externally applied voltage. Therefore the Fano factor 
7 can be computed as 



7 



2q\I\ 



(4) 



Before computing the ratio, the values of the numerator 
and of the denominator are uniformly averaged over the 
range of energy qV, assuming it much greater than kO 
(where k is the Boltzmann constant and 9 is the absolute 
temperature). In particular, our simulations have been 
performed using 500 values of energy in a range of 40 /leV 
around 9.03 meV. 

In Fig. [1] we show, with solid lines, the values of the 
Fano factor obtained for an 8 fim wide structure made 
up of a series of identical barriers. In particular, the 
reported results are relative to 0.425 nm wide barri- 
ers, with heights equal to: 0.8 eV (squares), 0.25 eV 
(triangles) and 0.07 eV (circles). These barriers (at 
the considered Fermi energy) have an average trans- 
parency r — J2n=i ksnP/-^ equal to about 0.1, 0.5 
and 0.9, respectively. In these simulations the dis- 
tance between adjacent barriers has been assumed to 
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FIG. 2: Same quantities as in Fig. [1] but averaged over 50 
sets of interbarrier distances. 



he D + 5, where D = 3 fim, and (from left to right) 
6 = 0, 10, -6, 3, -3, -9, 5, -10, -2 nm. 

We see that, as expected, for all transparencies the 
Fano factor for a single barrier is about 1 — F (not ex- 
actly 1 — r because different modes experience different 
transparencies and F is only an average value). For 2 
barriers our results are still in agreement with the semi- 
classical model of Ref. 7] and, specifically, with the re- 
sults of its Eq. (17). For more than 2 barriers we notice 
a sharp divergence from the semiclassical prediction and 
that no asymptotic 1/3 value is reached. Indeed, our re- 
sults show some dependence on the interbarrier spacings, 
but the overall behavior is already captured by the plots 
of Fig. [TJ A marked difference is observed only in the 
case of equidistant barriers, in which strong resonances 
between the different interbarrier regions play a major 
role, a case that we do not address in detail in this letter. 

We have also repeated our simulations making the 
same simplification adopted (for analytical convenience) 
in Ref. 0] , i.e. assuming a barrier transparency indepen- 
dent of the orthogonal wave vector of the impinging par- 
ticle. In detail, we have replaced the previously indicated 
scattering matrix for a barrier with that of an artificial 
barrier in which = — F and tb„ = V^, with 

F being the wave vector-independent transparency. No 
significant variation is observed in the Fano factor when 
such a change is included in our calculation (see dashed 
lines in Fig. [T]) . 

In order to remove the dependence of our results on 
the actual choice of the set of lengths of the interbar- 
rier regions, we have performed an average over several 
sets [13| . It has been shown that a similar approach is 
able to reproduce many effects of dephasing on trans- 
port [l3| . The results obtained by averaging over energy 
values and over 50 different sets of interbarrier distances 
are shown in Fig. [5] as a function of the number of cas- 
caded barriers, for the same transparency values as in 
Fig. [U assuming either a realistic barrier model (solid 
lines) or a wave vector-independent transparency (dashed 
lines). Also in this case there is no clear convergence to 
a common value of 1/3. 



3 



If we consider a situation with cascaded identical bar- 
riers characterized by transparencies that are indepen- 
dent of the wave vector, and average over random phases, 
all propagating modes give the same contribution to the 
noise behavior, because the different values of the longi- 
tudinal wave vectors are made ininfluent. Therefore it 
becomes possible to perform an analytical calculation of 
the Fano factor by considering a single mode and inte- 
grating over the phase of r„ between and 27r, for each 
interbarrier region. The analytical treatment can be car- 
ried out for the cases of 2 and 3 cascaded barriers (with 
transparency F), for which we obtain, respectively: 

2(1 -F) , 3(4-8F-K5F2 -F^) 
72 - 77^-FT7 and 73 = ^.t^ , nr^^ (5) 



(2-F)2 



16 - 24F + 9F2 



which are in agreement with the numerical results rep- 
resented with dashed lines in Fig. [2l Coherently with 
our previous discussion, the result for 2 barriers coincides 
with that from the semiclassical model of Ref. , while 
that for 3 barriers does not. This is consistent with the 
conclusions by Forster et al. [l^, who show that only 
in the case of a single probe (in our structure a probe 
should be included between each pair of barriers), cur- 
rent fluctuation statistics do not depend on the nature 
(phase averaging, elastic dephasing, or inelastic) of the 
probe itself. 

The key difference between the semiclassical and the 
quantum model consists in the fact that a semiclassical 
model (unless very peculiar assumptions are made [lB|) 
lacks localization as a result of complete incoherence, 
while a quantum model does exhibit strong localiza- 
tion [13] ■ In particular the 1-dimensional nature of the 
disorder represented by the randomly placed barriers 
makes the system effectively 1-dimensional, regardless of 
its actual dimensionality. In this case no mode mixing is 
introduced, and therefore the localization length is of the 
order of the mean free path; localization occurs beyond 
this length. Instead, in the case of 2- or 3-dimensional 
disorder, as in Ref. 15|], strong mode mixing makes the 
localization length approximately equal to the mean free 
path times the number of propagating modes [18]. 

In addition, in the absence of mode-mixing, it is not 
possible to consider the interbarrier regions as quasi- 
reservoirs, characterized by a well-defined occupancy, 
that depends only on the energy. This assumption is 
at the basis of the calculations of Refs. [3, d, [^, as well 
as of the semiclassical Monte Carlo numerical simulation 
by Liu et al. [l^ . 

Instead, in the absence of mode mixing, only a mode- 
dependent occupancy can be defined; we have computed 
it for the same structure as in Fig. [1] in which all pos- 
sible electron states can be divided into two sets: those 
injected from the left lead and those injected from the 
right lead. Therefore, if we define as ipnL and ipna the 
electron wave functions in the generic interbarrier region 
n resulting from an injection of the n-th mode (with unit 
probability current) from the left or right lead (respec- 
tively), the occupancy /„fj in the region for the n-th 
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FIG. 3; Values of the occupancy for 8 modes in the 5 interbar- 
rier regions of a series of six unequally spaced tunnel barriers 
with r — 0.1. Inset: distribution of the occupancy for the 
propagating modes in the 3rd interbarrier region. 



mode can be expressed as the ratio of the partial density 
of states related to injection from the left to the total 
density of states for that mode [20| : 



fnsi — 



In IV'"!. ]^ dxdy 



In \iPnL \ dxdy + ]i/'„^ ] dx dy 



(6) 



Results for the occupancy in the 5 interbarrier regions 
of a series of 6 unequally spaced barriers with an average 
transparency F — 0.1 are reported in Fig.[3]for a selection 
of 8 of the 320 propagating modes. It is apparent that 
these occupancies assume quite different values, with a 
strong dispersion that clearly appears in the inset, where 
we present the distribution of the occupancy in the re- 
gion between the 3rd and the 4th barrier. Therefore the 
assumption of quasi-reservoir behavior of the interbarrier 
regions is definitely not valid in this case. An exception 
is confirmed for the case of just two barriers (thus with 
a single interbarrier region), in which the occupancies 
are all equal and corresponding to the value predicted by 
semiclassical models. 

In the presence of mode-mixing, instead, the localiza- 
tion length Li is approximately equal, as already men- 
tioned, to the product of the elastic mean free path Lq 
by the number of propagating modes iV; therefore there 
can be a range of device length Ld values in which the 
condition for diffusive transport (Lq NLq) is 

satisfied, and thus the Fano factor can possibly reach the 
value 1/3 (as in the case of 2-d or 3-d disorder). 

An adjustable amount of mode-mixing can be intro- 
duced by applying a magnetic field orthogonal to the 
plane containing the device. We have computed the Fano 
factor for a 1 ^m wide structure with a series of 10 un- 
evenly spaced barriers, with average interbarrier distance 
500 nm. The results, obtained averaging over a number of 
interbarrier distance sets (40 for F — 0.1, 30 for F = 0.25, 
20 for F = 0.5, 10 for F = 0.75 and F = 0.9) are reported 
in the inset of Fig. [U for different choices of the bar- 
rier transparency F. In detail, the barriers are 66 meV 



4 



1.4 
1.2 

o 1 

S 0.8 



\ r=0.75 


r= 




^ r=o.9 


0.04 0.08 0.12 

Magnetic field (T),, Y= 




10 15 20 25 30 
Number of barriers 



35 40 45 



FIG. 4: Values of the Fano factor obtained for a series of 
identical barriers in a Ifim wide structure with an orthogonal 
magnetic field B = 0.1 T, represented as a function of the 
number of barriers. In the inset: Fano factor for a series of 10 
identical barriers, as a function of magnetic field. The dotted 
lines indicate the diffusive limit of 1/3. 



high, with a thickness of 4, 2.6, 1.56, 0.85 and 0.4 nm 
for the five considered transparencies; the Fermi energy 
is 9.03 meV. We observe that, after a quick drop, as the 
magnetic field increases the Fano factor settles around 
values that depend on barrier transparency. In the main 
panel of Fig. U] we report the Fano factor for a constant 
magnetic field of 0.1 T as a function of the number of 
barriers. We notice that for large values of the trans- 
parency, and thus large values of the mean free path, 
a diffusive transport regime (with a Fano factor of 1/3) 
is achieved only for a length much larger than the mean 
free path. On the other hand, for low barrier transparen- 
cies, and therefore reduced mean free path, localization 
effects appear before reaching the diffusive limit. To pre- 



vent this, we should significantly increase the number of 
propagating modes and thus the localization length. 

The issue is then whether in a practical system con- 
taining cascaded barriers large enough mode mixing takes 
place. Besides magnetic field, possible mechanisms lead- 
ing to mode mixing are scattering with irregularities in 
the potential (2-d or 3-d disorder) or phonon scattering. 
Scattering due to a disordered potential landscape can 
well lead to full mode mixing, but in such a case the 
Fano factor of 1/3 characteristic of diffusive transport is 
achieved anyway, independent of the presence of the bar- 
riers and cannot therefore be specifically attributed to 
their action. As far as phonon scattering is concerned, it 
can in principle introduce mode mixing, but in the pres- 
ence of strong phonon interaction the transport regime 
would not be the one we are interested in, and thermal 
noise would prevail. 

A Fano factor of 1/3 might in principle also be recov- 
ered, irrespective of the degree of mode-mixing, in the 
presence of a hypothetical elastic mechanism capable of 
suppressing phase coherence completely. 

The relatively large values of the phase coherence 
length that can be achieved in modern materials at low 
temperature and the low-field results for the Fano fac- 
tor presented in Ref. [l3| lead us to the conclusion that 
a superlattice or a series of electrostatically defined bar- 
riers in a channel containing a high-mobility 2DEG are 
more likely to exhibit localization and a Fano factor as 
predicted by our model, rather than a diffusive behavior. 
Numerical approaches along the lines we have presented 
could be instrumental in designing further experimental 
tests, which should be performed on structures with un- 
evenly spaced barriers. 

We are indebted to Prof. C. W. J. Beenakker for useful 
discussion. 
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